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1 Abstract 


Fast entangling gates for trapped ion pairs offer vastly improved gate operation times relative to 
implemented gates, as well as approaches to trap scaling. Gates on a neighbouring ion pair only 
involve local ions when performed sufficiently fast, and we find that even a fast gate between 
a pair of distant ions with few degrees of freedom restores all the motional modes given more 
stringent gate speed conditions. We compare pulsed fast gate schemes, defined by a timescale 
faster than the trap period, and find that our proposed scheme has less stringent requirements 
on laser repetition rate for achieving arbitrary gate time targets and infidelities well below 10“^. 
By extending gate schemes to ion crystals, we explore the effect of ion number on gate fidelity 
for coupling two neighbouring ions in large crystals. Inter-ion distance determines the gate time, 
and a factor of five increase in repetition rate, or correspondingly the laser power, reduces the 
infidelity by almost two orders of magnitude. We also apply our fast gate scheme to entangle 
the first and last ions in a crystal. As the number of ions in the crystal increases, significant 
increases in the laser power are required to provide the short gate times corresponding to fidelity 
above 0.99. 

2 Introduction 

Trapped ion systems have performed quantum algorithms on small scales CHS], however achiev¬ 
ing results inaccessible to the classical regime remains a significant challenge. There are vari¬ 
ous entanglement protocols [BHIl], many of which have been implemented CSH3I], where the 
required interaction between ions is mediated through the motion of the ions. The quantised 
motional modes can be used as an information bus, where spin-state information can be mapped 
onto the motional state, and subsequently onto another ion’s spin-state. This idea was outlined 
in the original proposal for quantum information processing (QIP) using trapped ions [6]. Gates 
mapping information between spin and motional states must be performed in the regime where 
spectral resolution of sidebands can be performed. This restricts the coupling strength to much 
less than the trap oscillation frequency, fl <C [32]. 

QIP schemes that offer advantages over their classical counterparts require large ion and 
operation numbers: at least tens of qubits and thousands of operations for quantum simulations 
|33| . The growing complexity of the motional spectrum with the number of trapped ions is a 
scaling problem for resolved-sideband gates. The timescale of entangling gates increases relative 
to the trap frequency and decoherence timescales, which are on the order of seconds or minutes 
depending on the internal states [341135] . The weak-coupling regime, where gate times with 
just two trapped ions are limited to much longer than the trap period, on the order of l//s, is 
thus prohibitive for large-scale QIP operations in a single trap. Scaling mechanisms using ion 
shuttling between traps and 2D architectures have been proposed, and promising steps taken in 
their implementation |10ll36Iffl!] . We explore the fast gate contribution as a scalable entangling 
operation within a single trap, which can be usefully combined with other scaling methods. 

Fast gates were proposed by Garcfa-Ripoll, Zoller and Girac in 2003 |l2] as a mechanical 
rather than spectral method for entangling ions. Here a controlled-phase gate is performed on 
two ions using state-dependent forces. Instead of resolving particular motional sidebands, fast 
gates excite multiple motional sidebands in the strong coupling regime where Q ^ v, offering 
gate times no longer limited by the trap frequency. The key element of an ideal fast gate is the 
motional control of the ions. State-dependent forces drive the motion such that the spin-states 
acquire relative phase and are maximally entangled. 

By moving into a new gate time regime, fast gates open new opportunities for scaling |43I149j . 
The application of fast gates to two ions in a multiple-ion crystal was first outlined by Duan [43] . 
Schemes faster than the local ion oscillation frequency are shown to entangle the ions, with 
arbitrary fidelity through multiple ‘cycles’ of the fast gate scheme. Instead of exciting collective 
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modes of an ion crystal, which would leave electronic and motional states entangled at the end 
of a fast gate addressing just two qubits, [l5] shows that only local modes are involved in the 
gate. More precisely, gates exceeding the local ion oscillation frequency only need consider the 
motion of neighbouring ions to those entangled by the gate. This locality allows a small number 
of degrees of freedom to determine a high-fidelity two-qubit gate scheme addressing two ions in 
an arbitrary length ion crystal. We explore this principle using pulsed fast gates, and extend 
fast gates to non-neighbouring ions. 

In section [3] we present the fast gate mechanism and conditions for the performance of a two- 
qubit gate, as well as the proposed pulsed gate schemes. Recent progress has been made towards 
fast gates using high repetition-rate pulsed lasers [T7lH8l 15(111^ . including implementation of 
the pulse pairs required for the gate [511152] and an exploration of spin-motion entanglement in 
the strong coupling regime [48] . This highlights the need for the optimisation of implementable 
pulsed fast gate schemes, which we provide in section 01 Here we compare the gate times and 
fidelities of the schemes given particular repetition rates, and provide the optimal regimes for 
each method. Applying this knowledge of fast gates to ion crystals is the main thrust of this 
paper; in section Owe explore gate hdelities and laser requirements for entangling neighbouring 
ions in a long ion crystal. We also present our results that a distant, non-neighbouring ion pair 
can be coupled in the fast regime using large momentum transfers, with motional restoration 
from simple symmetries in the pulse scheme. 


3 Fast gate formalism 

In this section we present a review of published fast gate schemes. 


3.1 Ideal gate and two-qubit conditions 

Fast gate schemes proposed with laser pulses use pairs of counter-propagating vr-pulses to give 
state-dependent momentum kicks to the ions, exciting various motional modes. The kicks 
are performed fast with respect to the ion trap dynamics. While detuned pulses can be used 
to drive phase trajectories using Stark shifts as in IISl, we focus on resonant transitions for 
simplicity. These state-dependent trajectories correspond to state-dependent paths in phase 
space that provide a particular phase to the given state. Area enclosed by a phase-space 
trajectory corresponds to the acquired phase m- We will consider particular phase-space 
trajectories, where we map the displacement and rotation of the centre of a coherent state. 
When the relative phase is vr/4, we have performed the desired gate operation according to the 
ideal unitary: 


f/j = i^pTca^ap 

Here is the Pauli Z-operator on the fcth ion, Vp and Op are the frequency and annihilation 
operator for the pth motional mode, and Tq is the total gate time. Ions one and two are targeted 
by this gate operation, however in general the two target ions can be chosen as desired. 

The product of rotation operators corresponds to the free evolution of the motional modes 
during the gate operation, over time Tq. There are L modes corresponding to the number of 
ions in the trap. At the end of the ideal gate, the motional modes are restored, meaning that 
they have been rotated by the free evolution that would have occurred without a gate. This 
restoration means that no heating is introduced by the operation. 

As presented in [121117], the two-qubit {L = 2) conditions for acquiring the desired phase 
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and restoring the motional conditions are given by 
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where N gronps of pulse pairs are applied with the group of pairs containing pulses 
arriving at time tk, and 6tmk = tm — tk > 0. The sign of Zk is the direction of the first pulse 
in each pair, and the trap centre of mass frequency and Lamb-Dicke parameter are v and rj 
respectively. Setting Cc and Cr to 0 represents motional restoration and 0 = vr/d represents the 
relative phase condition. More general conditions for arbitrary numbers of ions in the crystal 
are derived in section 15.11 in detail. 


3.2 Gate schemes 

Now that the conditions for an ideal fast gate operation have been established, we review 
proposed pulse schemes for applying the gate to a pair of ions. A general optimization search 
for the optimal pulse scheme given laser parameters and complete pulse timing and direction 
freedoms is a prohibitively difficult problem. Any restriction on the degrees of freedom for 
a more tractable search should ideally provide gate schemes with both an effective scaling of 
gate time with applied momentum, and a robust solution even for large numbers of ions. The 
following schemes use the motional symmetries discussed in Appendix 1 for motional stability 
and a tractable search space. We will show that each scheme has the optimal scaling of gate 
time with the momentum applied to the ions. 

3.2.1 GZC scheme 

The scheme proposed by Garcfa-Ripoll, Zoller and Cirac (GZC scheme) [42114^ is characterised 
by instantaneous groups of pulse pairs z sent at times t, interspersed with free evolution: 

z = (—2n, 3n, —2n, 2n, —3n, 2n) 

t = (-n, -T2, -T3, Ts, T2, Ti). 

At time —ti, 2n counter-propagating pulse pairs are applied along the trap axis (aligned 

with the z axis) to provide a Anhk momentum kick in the —z direction. The integer n determines 

— 2/3 

the gate time Tq, which scales with the total number of pulses in the scheme Np as Tq oc Np ' . 
For two ions this scheme exactly solves the condition equations ([2]),([3]) and (HD- 

The unitary kick at a given time can be written as the product of unitary kicks on each 
of the modes. In Figure [T] we show the scheme applied to two ions with n = 10 and a total 
gate time of 222 ns, in rotating and non-rotating frames for both motional modes. In rotating 
phase space, the initial motional state is stationary under no external evolution and the kick 
direction depending on the time of the kick and the mode frequency. Each 

mode-specific unitary has a kick strength that depends on the mode and the initial internal 
state of the ions. The internal states determine the phase space trajectories: two ions sharing 
the same state are displaced in the centre of mass frame and invariant in the stretch mode, 
while the reverse is true for two ions with different internal states. 

Figure [T] also shows the effect of laser limitations, where instead of instantaneous pulses, 
each pulse pair is separated by the laser repetition period. Instantaneous pulses at time ri, for 
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instance, become a gronp of 2n individual pulse pairs separated by the repetition period Tr and 
centred in time at ri. For large numbers of pulses, or very short gate times, a limiting factor for 
the gate is avoiding temporal overlap for the different gronps of pulse pairs separated by the finite 
repetition period. As the laser repetition rate slows, the approximation to the instantaneous 
ideal pulses fails and the gate fidelity can drop significantly. In this paper we focus on high 
repetition rates npwards of 300 MHz [50], which we will show still provide remarkably high 
hdelities. 


(a) 


COM 



Position 


(C) COM 
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(b) SM 



Position 


(d) SM 



Rotating position 


Rotating position 


- Ideal - fr 5GHz 


Figure 1: GZC scheme in non-rotating (a,b) and rotating (c,d) phase space with n = 10, for 
the centre of mass (a,c) and stretch (b,d) modes. Trajectories for the centre of a coherent state 
are plotted nsing dimensionless position and momentum. We show the performance of the ideal 
scheme, with an implied infinite repetition rate, and the performance of the scheme with a laser 
repetition rate fr = 5GHz. For the centre of mass plots, both ions are in the excited internal 
state. For the stretch mode plots, one ion is excited and the other is in the ground state. 


3.2.2 Fast Robust Antisymmetric Gate scheme 

Using the antisymmetric form of the GZG gate scheme, a search for the optimal gate scheme 
for implementing fast gates by splitting pnlses provided the following scheme m, which we call 
the Fast Robust Antisymmetric Gate (FRAG): 

z = (—n, 2n, —2n, 2n, —2n, n) 

t = (-n, -T2, -rs, Ts, T2, Ti). 

Like the GZC scheme, it is an exact solution to the condition equations (|2|), (|3|) and m 
for two trapped ions, nsing the timing freedoms for motional stability. Fignre[2]shows that the 
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phase-space trajectories closely resemble those of the GZC scheme, and again the effect of a 
5 GHz laser repetition rate is shown. This scheme was designed for splitting large-area pulses, 
but we will explore other advantages of this scheme for vr-pulse pairs directly applied to the 
ions. 




Position Rotating position 

- Ideal - fr5GHz 


Figure 2: FRAG scheme in non-rotating (a) and rotating (b) phase space with n = 10, for the 
centre of mass mode with dimensionless position and momentum. The ideal scheme and one 
with a laser repetition rate of fr = 5GHz are compared. Both ions are in the excited state. 


3.2.3 Duan scheme 

A simplified version of the scheme presented in [33], is described by 

z = (n, -2n, n) 
t = (0, Ti, 2 ti), 

where we have assimilated pulses with the same direction into instantaneous kicks to match the 
form of GZC and FRAG. As presented in [33|, the scheme was designed for pulsed lasers with 
finite repetition period Tr such that each pulse pair is separated by this period: 

z =(1, 1, 1, ..., 1, -1, ..., -1, 1, ..., 1) 

t = (0, Tr, “^Tr, ■■■, Ti Tr, Ti, ..., 3ti Tr, 3ti, ..., 4ti). 

The scheme was designed for minimal gate time, such that the repetition period between pulses 
is the only free evolution, and the scheme forms a triangle-like shape in rotating phase-space, 
as shown in Figure [3| 

This scheme takes the simplest form of a symmetry discussed in Appendix 1 to complete the 
triangular loop, using a single phase-space loop for maximal area rather than multiple triangular 
shapes as in the star-shaped schemes we have seen. Doubling the scheme, so that there are two 
phase-space cycles, cancels the first order of motional error. Here the scheme is simply repeated 
such that every pulse pair gives the opposite direction momentum kick to those in the first 
cycle. Figure [3] shows the doubled scheme, which takes almost twice as long as the single cycle 
for the 5 GHz repetition rate. 

Even with multiple cycles, the Duan scheme involves only one free variable in time, which 
is chosen such that the phase equation is satisfied. The scheme is thus not an exact solution 
to the condition equations. Accordingly, the symmetries employed provide gate speed require¬ 
ments such that the error terms are suitably small for motional stability. Multiple loops or 
cycles, appropriately designed, increase the motional stability as found in Appendix 1. We will 
construct the scheme with different numbers of cycles to enclose the approximate desired phase, 
then optimise to explore the tradeoff between fast phase acquisition and fidelity. 
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Figure 3: Duan scheme in non-rotating (a,c) and rotating (b,d) phase space with n = 10, for a 
single-cycle (a,b) and double-cycle (c,d) scheme. The centre of mass (COM) and stretch mode 
(SM) trajectories are shown. For the COM plots, the ions are in their excited states, while 
for the SM plots one ion is excited and one in the ground state. The laser repetition rate 
is fr = 5GHz, which provides the phase space curvature of the scheme. For the single-cycle 
scheme, ri = 60rr, while for the double-cycle scheme ri = 5lTr such that the phase gate is 
performed. Both schemes are not quite restored to their initial motional state, most clearly 
seen in (b) and (d), however the double-cycle scheme in (d) is more robust. 

4 Two trapped ions: gate optimisation and application 

The three presented schemes provide solutions to the fast gate condition equations ([2]), ([3]) and 
([^ for two ions in a trap. We use simulated annealing to optimise the pulse timings and the 
number of pulses n to maximise the gate fidelity and minimise the gate duration for a given 
repetition rate. We compare the resulting gate times and fidelities for each scheme. 

4.1 Fidelity measure 

Our performance measure is the state-averaged fidelity, derived in Appendix 2, which depends 
on the number of ions. For two ions, 

Faave = (e + cos(<))')) , (5) 

where x and y are the displacements in phase space for each mode, cj)' is the actual relative 
phase difference minus the desired tx/2 for different initial states. There is also dependence on 
the initial motional state, since 

( 6 ) 
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for the mean mode occupation n, which is set to be equal for each mode. We will vary the 
initial mode occupation n, and explore the effect of trap temperature on fidelity, in section 321 
The fidelity measure takes values between 0.5 and 1 due to the state-averaging effect, for 
relative phase between 0 and vr. For different numbers of ions, the fidelity function becomes 
more complex. However, it has the same exponential dependence on the initial motional state 
and the final phase space trajectory displacement from the ideal, as well as the sinusoidal phase 
dependence. 

4.2 Scheme performance for 

We use the ^^Ca'*' ion for our fast gate analysis as a typical ion candidate. Our parameters are 


u = 2 tt X 1.2 MHz 

( 7 ) 

A = 393 nm 

(8) 

T] = 0.16 

(9) 

re-i = re 2 = 0.1 phonons. 

(10) 


where we use the 5'i/2 to -P 3/2 transition to give the ions momentum kicks. The computational 
excited state, the metastable £> 5/2 state, is untouched by the resulting momentum kicks. Since 
the forces are only applied to the computational ground S' 1/2 state, this differs from the above 
derivation where equal and opposite kicks are assumed for the two states. The adjustment 
results in an effective halving of the momentum kick size in the condition equations derived, as 
shown in Appendix 3. 

For the given ion parameters. Figure 3] shows the gate time as a function of laser repetition 
rate, defined here as the rate at which a pair of counter-propagating 7r-pulses can be applied 
to the ions. Higher repetition rates correspond to stronger forces displacing the ions from their 
equilibrium positions, and the relative phase is acquired more quickly. Low pulse numbers in 
a scheme correspond to slower gates, such as the re = 5 (70 pulse pair) GZC scheme and the 
re = 9 (90 pulse pair) FRAG scheme for a 300 MHz repetition rate. Each scheme in the Figure 
has the optimal scaling of the repetition rate Tq oc fr as the repetition rate increases. This 
scaling is from the optimal gate time with total pulse number relationship [44] 

Tg oc iV-2/3, (11) 


with the approximation that Np = Tcfr- 

The high-fidelity schemes (FRAG and GZG) and the Duan schemes are marked by solid lines 
and dashed lines respectively in Figure 31 The high-fidelity schemes have error (1—Fidelity), or 
infidelity, below 10“®, while the Duan schemes have much lower fidelity as shown in Figure [5l 
The number of cycles of the Duan scheme is marked by Roman numerals in Figures 31 and 31 

Faster repetition rates lead to higher gate fidelities, as shown in Figure 31 The faster gate 
times simplify the motional conditions, as discussed in section 15.21 More Duan cycles also 
increase the robustness of the scheme to motional error. However, the gate time for four or 
more cycles of the Duan scheme was shown to be longer than the FRAG scheme. 

High motional occupation enhances the error according to the exponential terms in the 
fidelity function, up to the maximal error asymptote at 0.5. For very high fidelity gates, such 
as the two-ion GZG scheme with finite repetition rates, the infidelities for mode occupations of 
1 and 10 were found to remain on the order of 10“®. 

We will focus on the GZG and FRAG gates, with very high fidelities, as we extend the fast 
gate schemes to target two ions in ion crystals of varying length. The Duan schemes with one or 
two cycles may be faster, but have insufficient fidelity for QIP protocols with foreseeable laser 
repetition rates. 
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• GZC 

■ FRAG 

♦ Duan I 
▲ Duan II 
T Duan IV 
O Duan VIII 


Figure 4: Gate time as a function of repetition rate for the presented schemes. The fits show the 
optimal scaling of gate time with pulse number. Solid lines are fit to the high-fidelity FRAG and 
GZG gates with error (1—Fidelity) below 10“®. Dashed lines are fit to the Duan schemes, which 
have much lower fidelity. The Roman numerals enumerate the number of cycles, or triangles, 
in the Duan scheme. 


5 Scaling ion number: gate optimisation and application 

Here we consider longer ion crystals, L > 2, and the efficacy of applying entangling gates to 
both a neighbouring and distant ion pair. We thus explore the impact of fast gates on the 
scaling problem for trapped ions. 


5.1 General gate conditions 

We can extend the two-ion gate conditions to more general conditions for ion crystals with 
arbitrary length. Fast gate evolution is composed of momentum kicks and the free evolution of 
the ion motion. The momentum kicks are assumed to be fast relative to the ion motion, and 
evolve the state according to 

Ukick = ( 12 ) 


for a pair of resonant 7r-pulses addressing ions 1 and 2, chosen arbitrarily. Here z is the number 
of counter-propagating vr-pulse pairs comprising the momentum kick, k is the laser wavenumber, 
xi and X 2 are the ion positions. We consider pulse lengths much longer than atomic transition 
periods ~ 10“^^ s to satisfy the rotating wave approximation. 

The ion position can be written in terms of mode displacements and ion-mode coupling 
coefficients [53] . 


Xi = 




(13) 


where the dimensional mode position operator Qp is given by the mode annihilation and creation 
operators: 

Qp = y(ap + 4). (14) 

Using the Lamb-Dicke parameter. 


Vp = 



(15) 
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Scheme 
—^ Duan I 

Duan II 
Duan IV 


Figure 5: Gate error (1-Fidelity) as a function of repetition rate for the Duan cycle schemes. 
The roman numerals mark the number of cycles (triangles) in the plotted scheme. More cycles 
reduce the error for a given repetition rate, and increasing the repetition rate increases the 
motional robustness of the scheme. The solid line is for mean initial motional mode occupation 
hi = h 2 = 0.1, while the dashed and dotted lines represent mode occupations of 1 and 10 
respectively. 


the momentum kick is composed of displacement operators for each motional mode p: 


= Ilp=iDpi-2iz{b'f>af + bf’a^)r]p). (17) 


Ukick = 0'2)’?p(“P+4) 

I h(p) 


p= 

tL 


where there are L modes for L ions in the chain, and the displacement operator, 

Dp{a) = exp[aaj, — a*ap], (18) 


commutes for each mode p. 

The free motional ion evolution corresponds to a phase space rotation for each mode: 

TJ — £3 iUpStj^apap 

^p,mot — ^ 


(19) 


where 5tk is the time between the fcth and (fc + l)th momentum kicks. 

The total gate evolution is described by momentum kicks with direction Zk interspersed with 
free evolution until the next momentum kick; 


Ugate = 

Cpk = 2zk{b^f^crf + 4^V|)r/p. 


( 20 ) 

( 21 ) 


where N is the total number of pulse pairs. 

For a given mode p, the products of displacements and rotations give a total unitary 

Up\a) = e*^p|Q;), 

acting on an initial coherent state |a). Here 

N 

a = — i ^ 

k=l 

N m—1 

^ ^ ^ ^ [^pm^pk R'® 

m=2 k=l 


N 

aY.Cpke-^^-^'^ 

. k=l 


( 22 ) 


(23) 

(24) 
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The internal state is left invariant. 

The state |a) has been displaced in motional phase space, then rotated according to the ideal 
unitary in equation ([T]). The phase acquired is caused by the ion excursion from equilibrium 
and subsequent different potential experienced, giving rise to phase evolution terms. 

The total displacement is by an internal state-dependent amount 


N 


Cp — f ^ ^ Cpk& 




(25) 


k=l 


which we desire to be zero for the ideal gate process. This displacement also provides the 
unwanted phase term, 


St) ^ 5 


(26) 


which does not provide relative phase conditional on the states of two ions, as in ([T]). It is state 
dependent, and thus gives a relative phase term for each mode. When the unwanted phase-space 
displacement is zero, this relative phase is also zero. For nonzero displacements, the single-ion 
phase terms can be cancelled with single qubit rotations if significant. 

The condition for ideal motional evolution for each mode is thus 


N 


0 — i ^ ^ Cpi^ e 


iuptk 


k=l 


N 


(27) 

(28) 


k=l 


We expand the remaining phase term for mode p, neglecting global phase terms without 
internal state dependence: 


N m—1 


Cp ~ ^ ^ ^ ^ CpmCpk SYn.{yp{tra ^fc)) 
m=2 k=l 

N m—1 

= ^ ^ ZynZksii\{vp{t^ - 4 )). 


m=2 fe=l 

The phase for L ions is given by the product of the phase from each mode, 

thus to obtain the ideal unitary in equation ([T]), we need 

L N m—1 

7 = ^ ^ ZmZk sin{up{tjn - tk)). 


(29) 

(30) 

(31) 

(32) 


p=i 


m=2 k=l 


Equations (j28p and (|32l) provide L -|- 1 conditions to perform a fast gate with L ions. Note 
that for two ions, our derivation follows that of 


5.2 Gate scaling 

Entangling gates that require resolution of particular motional sidebands become increasingly 
challenging as the mode density increases, which occurs as the number of ions in the trap goes 
up. Since the gate time must be much longer than the frequency splitting [5l], this causes these 
gates to slow dramatically. In contrast, fast gates do not couple to particular sidebands, however 
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the number of motional conditions for the gate are proportional to the number of modes, or 
ions, in the trap. Solving these equations independently becomes prohibitively complicated for 
a practical gate, which would ideally be independent of the number of ions. 

The fast gate timescale provides a solution to this problem. A gate for two neighbouring ions 
that is much faster than the local oscillation frequency of other ions in the crystal only needs 
to satisfy the condition equations for local ions. This idea is applied in [35], using amplitude- 
controlled segments from a continuous-wave laser to demonstrate the reduction in degrees of 
freedom required for sufficiently fast gates on neighbouring ions. In fact, as we demonstrate 
in Appendix 1, only two equations are required for suitably stable and fast pulse schemes, and 
motional symmetries of different forms protect against different orders of the motional error. 
A given gate scheme is motionally stable below a gate time that depends on the degrees of 
freedom and symmetries of the scheme. This scheme-dependent stability is evident for two ions 
in section [31 Faster gate times push each scheme towards motional stability. Faster gates also 
accomodate more ions in the crystal, with higher local oscillation frequencies of the ions. 

Performing gates on non-neighbouring ions, in contrast to the locality exploited above, will 
necessarily involve the motion of the intermediary ions. However, the same reduced motional 
conditions of Appendix 1 restore the motion of each ion, thus there is a regime where the gate 
can be performed robustly. In deriving the local oscillation frequency, each ion is assumed to be 
at equilibrium, and here we require momentum kicks sufficiently large to break this assumption 
and couple ions within a timescale shorter than the oscillation frequency. 

Gates coupling non-neighbouring or distant ions in a chain via optimal control of laser pa¬ 
rameters have been proposed, following the fast gate formalism. In |46ll55j . the ion coupling 
is mediated by transverse phonon modes using small numbers of continuous pulse segments, 
resulting in gates much longer than the trap period. We will explore the pulsed laser require¬ 
ments to perform gates between non-neighbouring ions with high hdelity and speed, as well as 
the impact on the intermediate ions. 

Two ions can be addressed in long crystals using lasers that couple the ions through trans¬ 
verse modes, and it is straightforward to adapt our formalism to this case. We consider instead 
the axial modes, assuming that ions not targeted by the gate can have the addressed internal 
states shelved in non-interacting states to allow the lasers to shine on or close to the crystal 
axis. 

Within a single trap, different gate schemes are required depending on the chosen ion pair. 
Coulomb repulsion leads to tighter grouping of the ions in the centre of the trap, shown in 
Figure (6] which leads to varying inter-ion coupling strengths between neighbouring ion pairs. 
Alternatively, an anharmonic trap such as in [SS] has constant inter-ion distances, which would 
prove useful for coupling different ions with the same gate scheme instead of adapting for varying 
distance. 

It is important to note that even with perfectly scalable entangling gates, a single trap 
cannot be scaled to include an arbitrary number of ions. A strong radial frequency confinement is 
required relative to the axial frequency confinement so that the ions are prevented from buckling 
to a zig-zag formation, as shown in Figure [6j Lowering the axial frequency and addressing the 
radial modes has been used to manipulate large numbers of ions in a single trap [56|. The 
anharmonic trap in |55| provides stable confinement for many ions, and fast gates have also 
been proposed for a 2D architecture using controllable continuous pulse segments [57j . While 
our pulsed gate analysis could be extended to this architecture, we focus here on standard linear 
traps. 

5.3 Coupling neighbouring ions 

We are interested in the scaling performance of pulsed gate schemes, as well as the laser re¬ 
quirements for particular benchmarks. We explore the performance of fast gates applied to the 
first two ions of a chain, as well as the effect of performing the gates on the central ion pair. 
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Figure 6: Equilibrium ion positions are shown in dimensionless units for different crystal lengths. 
For linear confinement of L ions, the ratio between radial and axial COM frequencies corlu^z > 
0.63L°'®®® [58l[59], and we plot the lower bound equality with the crystal length. 


'When ions are added to a trap with frequency i', with the frequency fixed and independent 
of the number of ions, the inter-ion spacing decreases. Increasing the number of ions in a trap 
thus increases the coupling strength of the target ions, and the gate is performed faster as shown 
in Figure [7l|a). 

The inhdelity grows with the number of ions in a trap with fixed frequency, as seen in Figure 
[7[c). Although added ions reduce the gate time, the motional requirements on the gate time 
for high fidelity also become more stringent. Here the number of pulses in the applied gate is 
chosen to provide the maximal fidelity, given the repetition rate and number of ions. Faster 
repetition rates give faster gates with higher fidelity, using more pulses. 

As shown in Figures 7(a) and 7(c), the GZC and FRAG schemes respond similarly as ions are 
added to the trap. However, the FRAG scheme is both faster and achieves higher fidelities. For 
the following results, we consider just the FRAG scheme for simplicity due to these advantages. 

Varying the trap frequency changes the separation of ions in a trap. In Figure [T^b), we show 
the gate time in terms of the separation distance of two ions. Here we have a fixed number of 
pulses N, and to reach the necessary relative phase, the phase equation p2|l determines that 
the inverse cubic dependence of the acquired phase on the ion separation is counteracted by a 
linear increase in pulse timings. The gate time thus scales linearly with the inter-ion distance, 
as shown in the hgure. 

The distance scaling obscures the effect of adding ions to the crystal. As ions are added, we 
can relax the trap frequency to maintain a constant distance, and coupling strength, between 
two ions. An optimal fast gate applied to these two ions has pulse timings that vary minimally 
with the number of ions in the crystal; the gate is ion-number independent. In Figure [HI we 
use this method to apply an optimal fast gate to the first two ions in the crystal, with different 
numbers of total ions. Fidelity decreases with the crystal length as the motion of each ion is 
not completely restored. While the distance between ions one and two is fixed, the distance 
between ions two and three decreases as more ions are added to the crystal, as can be observed 
from the ion crystal equilibrium positions in Figure [6l The ions neighbouring the two target 
ions thus become increasingly coupled to the motion of the target ions as the crystal grows in 
length. 

We similarly fix the distance of the middle two ions in the crystal and apply the fast gate 
to these two ions, also shown in Figure [HKa). Here the fixed distance between the middle, 
target, ions is smaller than the distance from a middle ion to its outer neighbour. The fidelity 
increase in the figure is caused by this decrease in motional coupling to neighbouring ions. This 
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(a) 


• FRAG scheme ■ GZC scheme 


Figure 7: Inter-ion distance has a strong ef¬ 
fect on gate time and fidelity, (a) Gate time 
for fixed trap frequency u = 1.2 x 2 tt MHz 
as a function of ion number in the trap. The 
gate is performed between the first two ions 
in the crystal, with a 5GHz repetition rate. 
The fidelity is above 0.99 for all save the 10 
and 20 ion GZC data points, with hdelity 
0.98 and 0.96 respectively, (b) Gate time for 
FRAG scheme as a function of the distance 
between two trapped ions. The same num¬ 
ber of pulses are used for each gate. The 
gate time depends linearly on ion separation, 
(c) FRAG (solid line) and GZC (dashed line) 
scheme error as a function of ion number, for 
ions added to a trap with fixed axial confine¬ 
ment u. The schemes are applied with an 
optimised number of pulses for each repeti¬ 
tion rate and ion number. 



difference is more pronounced for even numbers of ions, such that there truly are two middle 
ions in the crystal. As the number of ions increases, the middle ions and their neighbours in the 
crystal become approximately equidistant, and the motional fidelity reaches an asymptote. For 
gates between middle or end ions in the crystal, the repetition rate dependence is clear; higher 
fidelities correspond to faster gate times. 

The effect of gate speed on motional stability is clearly seen in Figure [SKb). The driven 
displacement of the local ions is restored effectively at the end of the faster gates, as the motional 
restoration equations reach their stable regime. This same behaviour is seen when the gate is 
performed on the middle ions of the crystal rather than the first two ions. As shown in [33], for 
gate times faster than the ion recoil frequency, only local phonon modes are excited and only 
neighbouring ions affect the gate operation. In Figure [HKb), there are five ions in the crystal 
and the local oscillation timescale of the third ion is 280ns. The gate times for the system, 
corresponding to different repetition rates, pass below this timescale and show correspondingly 
stable results in (a) and (b) of the figure. We will consider the displacement amplitude and the 
harmonic approximation in section [6j 
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Ion Number Time (ns) 

Repetition rate (GHz) 


Figure 8 : The distance between addressed ions is fixed, to clarify the ion scaling behaviour, (a) 
Motional infidelity as a function of ion number given a repetition rate. The gate takes 494 ns, 
293 ns, 154 ns, and 89 ns for repetition rates of 0.3 GHz, 1 GHz, 5 GHz and 20 GHz respectively. 
Solid lines are for the gate performance on the first two ions in the crystal, dashed lines represent 
the gate performance on the middle ions, (b) Driven displacement from free evolution of ions 
three (solid) and four (dashed) over the gate operation with different repetition rates, for 5 ions 
in the crystal. The 0.3 GHz repetition rate is neglected for clarity. Increasing the repetition 
rate reduces the displacement of neighbouring ions during and after the operation. The local 
oscillation period of the third ion in the crystal T 3 , neighbouring the target ions, is around 280ns. 
For repetition rates of 1, 5 and 20 GHz respectively, the gate times are I.O 5 T 3 , O.SSTs and 0.32r3, 
reflecting the stability gained by operating faster than this timescale. Both addressed ions are 
in the excited state. 
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5.4 Arbitrary couplings 


Now that we have explored the laser regimes for high-fidelity, fast gates between neighbouring 
ions in a chain, we consider non-neighbouring ions. It is possible to perform scalable operations 
using only neighbouring-ion operations and SWAP gates to couple distant ions. However, 
it would certainly be simpler if a direct entanglement operation between distant ions were 
achievable. We consider the laser requirements of pulsed fast gates between two distant ions. 

First, consider the mechanism for coupling neighbouring ions in a large crystal. High fidelity 
gates involve fast gate times such that only local modes are involved in the gate operation. In 
this regime, there can be no communication between distant ions via shared modes. 

The phase equation (|32p quantifies the interaction strength and the gate time required for 
the requisite relative phase. For distant ions, the coupling coefficients decay with the cube of the 
separation. Either a longer entanglement timescale or larger momentum transfers are required 
to excite non-local modes. The motional condition equations (1281) don’t change depending on the 
ions addressed, meaning that sufficiently fast gates will restore every mode, where the required 
timescale depends on the scheme stability as discussed. 

We approach the problem by first finding the optimal solutions to the instantaneous mo¬ 
mentum kick case, with an infinite laser repetition rate. The effect of a non-local gate on 
the ion crystal can be seen in Figure [9l where a FRAG gate entangles ions one and five in a 
five-ion crystal. Intermediate ions are key to the coupling of the distant ions, and accordingly 
have state-dependent motion. The driven motion of ions one and five is much greater than the 
driven displacement of the intermediary ions. To couple the ions with fidelity > 0.99, very large 
instantaneous momentum kicks are applied; n = 400, thus the largest kick uses 800 vr-pulse 
pairs. Evolution between the momentum kicks appears linear, showing that the free evolution 
of the ions dominates. In Eigure[9|^c), the restoration of the ions’ motion is shown to be effective 
for each ion. 

Our analysis of the GZC and ERAG solutions in the two-ion and neighbouring-ion cases, 
which used finite laser repetition rates, revealed that the fidelity was preserved remarkably 
well with fast repetition rates, at quantum computing threshold levels. We now apply finite 
repetition rates for our distant-ion gates. 

Schemes with finite repetition rates have a maximum momentum transfer in time, unlike the 
instantaneous pulse scheme of EigureO This can be seen, in the form of a maximum curvature, 
for the driven displacement of the ions. Eigure fTOl bl shows this effect in a displacement plot 
of three ions, where ions one and three are coupled, for three different repetition rates. The 
schemes are optimal given the repetition rates and gate time ~ 140 ns, and have comparable 
fidelities close to 0.98. The slowest repetition rate, 5 GHz, requires the largest momentum 
transfer to reach the necessary phase-space area corresponding to the target relative phase, 
given its lower curvature limited by pulse separation time. 

Eigure fTOl al shows the fast gate fidelity as a function of gate time, coupling ions one and 
three in a three-ion crystal. The infinite repetition rate corresponding to instantaneous large 
momentum kicks gives fidelity arbitrarily close to 1 as the gate time decreases and the motional 
fidelity of the scheme improves. As the gate time increases, the motional instability grows and 
the fidelity drops away. Einite repetition rates have poor performance when the gate is not long 
enough to acquire the target relative phase, as the repetition rate limits the total number of 
pulses in a given gate time. The optimal gate performance with a finite repetition rate occurs 
where the gate is long enough to entangle the distant ions with the desired phase, and short 
enough that the gate is motionally stable. Higher repetition rates reach a higher optimal fidelity 
with faster gate times. Lower repetition rates spread the applied force over longer time. This 
causes the motional stability to improve with increasing gate time, and be lower overall. 

Eigure m^) shows the gate fidelity for coupling ions one and five in a five-ion crystal. The 
effect of repetition rate is stronger for more ions; a repetition rate of 10 GHz is not sufficient 
for high fidelities. Even with an infinite repetition rate, achieving high fidelity requires much 
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Figure 9: (a) Driven displacement from free 
evolution of the five ions in the trap during 
the FRAG controlled phase gate operation. 
Ions one and five are entangled by the opera¬ 
tion, with fidelity of 0.99 using large instan¬ 
taneous momentum transfers with n = 400. 
Here ions one and five are in the same inter¬ 
nal state (both ground or excited), and we 
consider equal and opposite momentum kicks 
for each state for clarity in this displacement 
figure, (b) Ions one and five have different in¬ 
ternal states for these dynamics, \eg) or \ge). 
(c) Ions one and five have the same internal 
states for these dynamics. The final displace¬ 
ment of the ions from their initial positions 
is much less than 1 nm. 


lower gate times than for gates coupling ions one and two or one and three. 

6 Limitations 

The fast timescales of the entangling gates we have explored are achieved by careful control 
of momentum transfers. Imperfect vr-pulses lead to incomplete population transfer, leading to 
unwanted state changes, incomplete phase-space trajectories and motional heating. The timing 
of the pulses is also important, and timing errors also affect the final motional state and phase 
acquisition. General characterisation of these errors is limited, as they depend on the scheme 
applied and the number of pulses. As explored in m. the pulse timing is remarkably stable 
up to systematic shifts of beyond lOps, and this stability is reflected in the robust nature of 
the GZG and FRAG schemes under non-instantaneous laser pulses. In contrast, the pulse area 
can have a dramatic effect; in m the authors found that a 1% systematic pulse energy error 
leads to a worst case fidelity estimate of almost 1% for just a four pulse pair scheme. The worst 
case fidelity contains combinatorially increasing prefactors of the error in pulse area with the 
number of pulses. 

Initial motional states have no effect on perfect two-qubit gates, however imperfections or 
more qubits in the crystal introduce initial state dependence. The fidelity measure, equation (|5|). 
decays exponentially as the product of the mean vibration mode occupation and the final 
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Figure 10: (a) Optimal fidelity for a 
given gate time and repetition rate 
using the FRAG scheme to cou¬ 
ple ions one and three in a three 
ion crystal. The three compara¬ 
ble data points marked by a rectan¬ 
gle are explored in (b). (b) Driven 
displacement from free evolution of 
the three ions in the trap during 
the controlled phase gate operation. 
The different lines correspond to 
the different repetition rates of (a), 
for the marked optimal schemes 
with gate time close to 240ns, and 
fidelity close to 0.98. Ions one and 
three are in the excited and ground 
states respectively, with equal and 
opposite momentum kicks applied 
to the states for clarity of displace¬ 
ments. (c) Optimal fidelities for en¬ 
tangling ions one and five of a five 
ion crystal. 
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vibrational displacement error. This effect is quantified in Figure O where the impact of mean 
motional occupation is shown for imperfect two-qubit gates. Note also that gates performed in 
sequence will lead to an increasing error as the small heating from the first gate amplifies the 
error in the second gate, and so on. 

The driven motion of the addressed ions increases for faster gates or more distant ions. 
The lengthscale of this motion is still much smaller than the laser wavelength for the data we 
have explored, meaning that addressibility is no issue. However, the harmonic approximation 
made for the motional modes breaks as the ion displacement becomes significant relative to the 
ion separation. A driven displacement of 100 nm is still much smaller than the ~ 5 fim ion 
separation, and indeed oscillation amplitudes on this scale only change the harmonic spacing of 
the relative motion frequency by 0.2% [58]. Entangling distant ions requires increasing maximal 
displacements as more intermediary ions are added, and the inter-ion separation also decreases. 
As the current limits on laser repetition rates are overcome, increasing the possible momentum 
transfers, anharmonicities will become a limiting factor. Typical laser powers required for 
can be estimated using the 12 nJ vr-pulse energy presented in |l8|: average laser powers 
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of around 12 W and 60 W correspond to repetition rates of 1 GHz and 5 GHz respectively. 


7 Conclusions 

We have considered a general formalism for scalable fast gates using pairs of vr-pulses. Proposed 
pulsed gate schemes were analysed for different laser repetition rates, and the FRAG scheme 
was found to provide optimal gate times and fidelities. The relationship between repetition rate, 
gate time and fidelity was extended to neighbouring and non-neighbouring ions in long crystals. 
Increasingly large momentum transfers, corresponding to fast repetition rates, are required as 
ions are further separated in the crystal both to preserve motional robustness and to provide 
coupling faster than the local oscillation frequencies of the ions. We present here the repetition 
rates required to achieve information processing benchmark fidelities and times for two-qubit 
gates in arbitrary-length ion crystals. 
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We present three appendices to consider in detail the effect of symmetry on the motional 
conditions, the details of deriving our fidelity equations, and the adapted equations for treating 
ions with momentum kicks applied to only one internal state (‘asymmetric’ kicks). 


A Symmetry and the motional conditions 


Equation (25) in the main text gives the general condition for restoring motional mode p: 


N 

k=l 

where the final displacement of the motional mode is given by 

N 

Cp = -2ir,pibf'^al + , 

k=l 

N 

= isp ^ Zke^''p ^’‘, for 
k=l 

Sp = -27?p(6^^Vf + 62 ^V|). 

Expanding the exponential as a Taylor series in Cp, 

Cp — isp y ^ Zfc (1 + ii'ptk H— - —I—-—h ■■■), 

k=l 


and for Cp = 0, we require that 


Re(Cp) 

0 

Im(Cp) 

0 


N 




,'^Zk{-iyptk + + ...) = 0 


k=l 


N 


y^^zkjtk 


'^p^k 


k=l 


N 


,'^Zkil - + ...) = 0, 


k=l 


N 


k=l 


For 1 Vpt'l for each tk, the motional conditions reduce to 


(33) 

(34) 

(35) 

(36) 

(37) 

(38) 

(39) 

(40) 

(41) 


N 

0 = "^Zkik, 

(42) 

k=l 


N 


0 ^ ^ Zk 1 

k=l 

(43) 


with no mode dependence. This regime and these conditions correspond to free evolution, before 
the trap evolution has any significant effect. The required regime may be difficult to reach since 
the dominant term in our phase condition contains the third order of the pulse times, indicating 
that phase accrual has strong dependence on the trap evolution time. Strong momentum kicks 
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are required to accommodate shorter gate times, as explored in the main text. Close to this 
ideal regime, error in the scheme {Cp 7 ^ 0 ) will be mode dependent, and the dominant error term 
is from the sum in equation (I4U|1 . proportional to for each kick k. We can broaden the 

mode-independent regime by imposing constraints on the scheme pulse directions and timings. 
This provides increasing motional robustness. 

A simple example is to impose a reflected antisymmetry, such that the momentum kicks z 
and their timings t obey 

Z — (fll) ®2) I Q-m— 1) •••) Q-l) 

t — ( Tlj "^2; An I Anj An—Ij •••) A)) 

This sets Im(Cp) to zero by virtue of its even powers of tk, while we still require 


Re(Cp) = 2sp ^ ak sin{upTk) = 0, (44) 

k=l 

0 = .... (45) 

k=l 

The highest order term in our new motional condition equation corresponds to equation 
(j42|] . while we have subsumed equation (I43p into the symmetry condition. Assuming mode 
independence means that the dominant error term is proportional to aki'pT^, which is on the 
order of the cube of the pulse times. The asymmetry decreases the error size from a quadratic 
power of pulse times above, providing motional stability. 

We can introduce a more complex symmetry to further expand our region of mode-independent 
total gate times: 

z = (oi, ..., ai, I —ai, ..., —oi, I —ai, ..., —o;, | a/, ..., oi) 

i=(-/-A, \-f + Ti, ...,-f + Ti, |/-ri, \f + Tl, ...,f + Ti), 

for / > Ti > T 2 > ... > A. Here Re(Cp) is zero, and we require 

i 

Im(C'p) = 2sp ^ Ofc cos(i/p(/ -k Tfc)) - Ofc cos{up{f - Tk)) = 0, (46) 

k=l 
I 

0 = 

k=l 

where if the first term dominates the sum, we have a mode-independent condition equivalent 
to equation (1421) . The largest error term in the sum is on the order of on the fourth 

order of the pulse times for / ~ A. Again, the added symmetry conditions increase the motional 
stability region for the scheme. 

Adding ions reduces the mode-independent regime, since the mode frequency plays a part 
in the terms we require to be negligible, and these frequencies increase with ions added. Extra 
degrees of freedom found in the GZC and FRAG schemes relative to the Duan single degree 
of freedom also provide extra motional stability. Note also that while motional conditions can 
be satisfied independent of ion numbers for sufficiently fast gates, the phase equation still has 
mode dependence. However, only local phonon modes can be excited as discussed in the main 
text. 


’^pfTk “ -if - Tk)‘^) + 


(47) 
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B Fidelity measure derivation 

The computational fidelity before state averaging is given by 


Fi = Trm[{lpo\U^dUre\lpo){'lpo\ ^ PmU^eUidl'lpo)], 


( 48 ) 


where C/re and C/id denote the real and ideal gate operations, including motional displacement 
for the real gate as performed by our gate schemes to follow. The motional trace Tr^, is taken 
along with the inner product of the ideal and real operations with respect to the internal ion 
states. The initial internal state of the ions is {ipo), and the initial motional state is given by 
the density operator pm- In the computational basis {\gg), \ge), \eg), jee)}, the ideal unitary is 
given by 


Ui 


/ 0 0 0 \ 
0 0 0 
0 0 0 
\ 0 0 0 / 


(49) 


up to global phase. 

The real operator, representing our approximation to the ideal unitary with some error, is 
given by 


C/re 


/ 

V 




99 


0 

0 

0 


0 0 0 \ 

6 **^ 9 = 0 0 

0 0 
0 0 e*'^'='=i)ee y 


(50) 


representing phase changes (l>ij as well as motional displacements bij specific to internal states, 
b j £ {H) c}. Note that there are no off-diagonal terms as the initial internal states are preserved 
through the real operation due to the assumed perfect vr-pulse pairs providing momentum kicks 
while perfectly restoring state population. 

Recall that the final displacement is given in equation (|33l) for each mode p: 


N 

= (51) 

k=l 

The total displacement operator is the tensor product of each mode displacement. The internal 
state dependence comes from Cpk'- 

Cpk = 2Zk{l)f^<Tl + 4^Vf)??p, 

N 

Cp = -2i{h^\l + h^^Dpp ^ 

k=l 

We have defined the motional operators 

blf = {iij2\bp{Cp)\hj2), (54) 

for i and j internal states of ions one and two. Each mode displacement is related; for both 
ions in the same or opposite internal states these relations are simple: 

O// = (£>1?)' 
oj? = (4?)* 


(52) 

(53) 
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(55) 

(56) 






due to the opposite direction of each displacement kick. The other relations are by a scalar 
displacement amount, such that the motional displacement Cp for internal states |ee) and \eg) 
are related by 


hip) I hip) 

Cpiee) = -^^Cpieg). (57) 

The tensor product composite, mode independent, displacement operators for each state are 
related in the same way; 


Dpg = n^=^D^) = Dl (58) 

4e = = Dig. (59) 

We derive the phase acquired for each state from equation (24) in the main text, the mode- 
dependent phase term; 


N m—1 N 

Cp — I ^ ^ CpmC-pk sin(^'p(tm tk)l 3fi[Q; 'y ^ Cpk(^ ^ (60) 

m=lk=l k=l 

for an initial motional coherent state a. This second component of the phase is zero when the 
motion is restored at the end of the gate. As discussed in the main text, when the displacement 
is nonzero and this term is significant, we can correct for it using single-ion rotations as it is not 
entangling. The total phase term, from equation (30) in the main text neglecting global phase, 
reduces to; 


L L N m—1 

^ ^ sin(i/p(tm - 4))- (61) 

p=l p=l m=l k=l 

The presence of both [/re and its Hermitian conjugate in the fidelity equation ensure that global 
phase does cancel. Internal states determine the value of the sum, and we have defined 

L 

(pij = {hj2\^Cp\kj2), (62) 

p=i 


for ij G {g,e]. 

Total phase acquired by particular internal states is related according to 

4*gg — 036 — 063 — 066; (63) 

as each mode-dependent term in the sum changes sign for different internal states of ions one 
and two. 

For a general initial internal state 

|0o) = am\gg) + am\ge) + aio\eg) aii|ee), (64) 

A= (0o|4i[i[/re|0o) (65) 

+ |aii|2e*(‘^--^/^)T)e6. (66) 

Using the cyclic nature of the trace, we have 

Fi = TTm[A^ Apm]- (67) 
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The trace becomes an expectation of pairs of motional displacement operators with respect 
to Pm- Since operators on different motional modes commute, we can group the mode-dependent 
components, for example 

= ( 68 ) 

For a single mode, products of displacement operators are determined by 

b{a)b{h) = + h) (69) 

= b{a + b) b = ca for a scalar c. (70) 

We have seen that the arguments of our displacement operators, Cp from equation (|54l) . are 
indeed related by scalars for the same mode and different internal states. 

If the initial motional state pm is separable with respect to the motional modes, then the 
expectation values can be taken individually for separate modes. We assume an initial thermal 
product state as a typical motional distribution, 

OO 

p(p) = (1 - (71) 

n=0 

Pm = p^^ ® p^^ ® p^^ ■: (72) 


represented in the number basis, with temperature T and k the Boltzmann constant. 
The expectation value for a single mode with a displacement z follows; 

OO 

{ 0 ( 3 )),„ = (1 - 

n=0 


= g-bP(l/2+rip)^ 


(73) 

(74) 


where rip is the average phonon population of the pth motional mode. 

Choosing a particular internal state has limited use as a figure of merit for a quantum gate, 
which takes inputs of different initial forms in practice, the fidelity of each being important. We 
consider the average fidelity, taking an even weighting over all initial states. This corresponds 
to integrating over Fi for the possible values of the coefficients; an integral over the three- 
dimensional surface (given four basis states) of a hypersphere of radius one. We now have the 
formalism in place to calculate the average fidelity, which depends on the number of ions as this 
determines the phase and displacement equations. For two ions, 

i"2ave = (6 + , (75) 

where = (1/2 -|- hj) is a function of the motional mode mean occupation level, and Cp is the 
final displacement for mode p defined in equation (l53l) using equal and opposite internal states 
for nonzero Ci and C 2 respectively. 

For more ions, there are more complicated internal state dependencies in the Cp final dis¬ 
placements, so we define 

N 

Cp = 2r?p^Zfce*"^*S (76) 

k=l 


which is independent of the internal state for a simpler fidelity definition. For three ions, the 
fidelity becomes 


Faave = —(6-|-_|_ g-5.33mi|ci|2-2m2|c2p-0.66m3|c3|2 

_j_ 2(g-l-33mi|cip-2.66m3|c3p _|_ g-1.33mi|cip-2m2|c2p-0.66m3|c3p^ 2 ,^^^) ( 77 ) 

The average fidelity for no relative phase and very large motional displacement is 0.5, a point 
to keep in mind as we consider the fidelity values for the schemes in the main text. 
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C Asymmetric momentum kick equations 


If the laser providing the momentum kicks addresses a transition from one computational state 
to an auxiliary state, then the other computational state is untouched by the ‘kick’. Without 
loss of generality, we choose the ground state to undergo the momentum kick, while the excited 
state is left invariant. The unitary kick operator, equation (12) in the main text, becomes 

(78) 

where a?|0)j = 0 and af|l)i = |l)i. Thus (af)^ = a\ and a\ = (Ji + ct?)/ 2 for identity operator 

I. 

The mode basis transformation and the displacement and rotation operator expansion take 
place as for symmetric kicks in section 4.1 of the main article, with the adapted motional 
displacements 


N 

Cp = + bfal) (79) 

k=l 

(80) 

for each mode. The phase acquired for a given scheme becomes 
L L N m—1 

+ 4 ^^ “ 2 )^ X] X] ZmZkSin{Up{tm - tk)) (81) 

p=l p=l m=2k=l 

L N m—1 

= 4^??p((6i^V«i + (4^V«2 + 2b^fh^^''ala2) ^ ^ ZmZk sm{i^p{tm - tk)), (82) 

p=l m=2 k=l 

which includes non-entangling state-dependent phase. This causes deviation from our target 
gate, which can be corrected by single-qubit gates. 

We thus have the asymmetric phase condition equation 

L N m—1 

^ ^ ^ ZraZk Sm{Up{tm - tk)), (83) 

p=l m=2 k=l 

where the right hand side is one quarter the size of the symmetric kick phase equation (32). 

A symmetric scheme that solves the condition equations will thus provide only vr/lG phase 
in equation (l83jl . or 7r/8 relative phase between symmetric (|00), |11)) and asymmetric (|10), 
|01)) states. The adapted displacements and phase can be directly substituted into the fidelity 
derivation. 
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